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ABSTRACT 

We report the detection of 21-cm absorption at z a ^ s 



= 0.656 towards 1622+238 (3C336). 



The line is very broad with a Full- Width Half Maximum (FWHM) of 235 km s~ , giving a 
velocity integrated optical depth of J rdv = 2.2 ± 0.2 km s _1 . The centroid of the line is 
offset from that of the known damped Lyman-a: absorption (DLA) system by 50 km s -1 , and 
if the Lyman-a: and 21-cm absorption are due to the same gas, we derive a spin temperature 
of T s ^ 60 K, which would be the lowest yet in a DLA. The wide profile, which is over four 
times wider than that of any other DLA, supports the hypothesis that the hydrogen absorption 
is occurring either in the disk of a large underluminous spiral or a group of dim unidentified 
galaxies, associated with the single object which has been optically identified at this redshift. 

Key words: quasars: absorption lines - cosmology: observations - cosmology: early Uni- 
verse - galaxies: ISM - galaxies: individual (3C 336) 



1 INTRODUCTION 

Redshifted absorption systems lying along the sight-lines to distant 
quasars are important probes of the early-to-present day Universe. 
In particular, damped Lyman-a absorption systems (DLAs), where 
JVhi ^ 2 x 10 20 cm -2 , are useful as they account for at least 80% 
of the neutral gas mass density in the Universe (Prochaska et al. 
2005). Since the Lyman-a transition occurs in the ultra-violet band, 
direct ground based observations of this transition are restricted to 
redshifts of z > 1.7, although observations of the H I spin-flip tran- 
sitions at A rcst = 21 cm can probe neutral hydrogen from z — 0, 
thereby providing a useful complement to the high redshift opti- 
cal data. Furthermore, 21-cm observations can be used to confirm 
and complement reported variations in the fundamental constants 
over the history of the Universe (Murphy et al. 2003). Unfortu- 
nately, redshifted radio absorption systems are currently very rare 
with only 50 Hi 21-cm absorption systems known (summarised 
in Curran et al. 2006), 17 of which are DLAs, with only nine of 
these having sufficiently high quality UV/optical and 21-cm data 
(Tzanavaris et al. 2007). We are therefore currently undertaking a 
survey of, as yet unsearched, DLAs occulting radio-loud quasars 
in order to increase the number of 21-cm absorbing systems and 
in this letter we report our first detection, in the intriguing DLA 
towards 1622+238 (3C336). 
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2 OBSERVATIONS AND DATA REDUCTION 

The observations were performed on 9-10 luly 2006 with the West- 
erbork Synthesis Radio Telescope (WSRT) as part of our survey 
(described in Curran et al. 2007). The front-end was the UHF-high 
receiver tuned to 857.68 MHz (or z a b s = 0.6561, Rao & Turn- 
shek 2000), backed with a band-width of 5 MHz over 2048 chan- 
nels (dual polarisation), giving a channel width of 0.85 kms -1 . 
The two orthogonal polarisations (XX & YY) were used in order 
to allow the removal of any polarisation-dependent radio frequency 
interference (RFI). After the flagging out of time-dependent RFI, 
12.0 hours of good data remained, although there was some RFI 
remaining on some baselines, particularly in the XX polarisation. 
After further flagging, 63 full and partial baseline pairs remained. 
Additionally, frequencies below 856.6 MHz (channels above 1500) 
were removed due to a steep 1 ly dip in the flux. 

The data were reduced using the MIRIAD interferometry re- 
duction package, with which we extracted a summed spectrum 
from the emission region of the continuum map. Upon this, we 
found absorption towards the peak of the radio continuum appar- 
ent in both polarisations. After subtracting the continuum with a 
first order polynomial from the line free channels of the combined 
polarisations, a single Gaussian gives the fit shown in Fig. 1. 

The obvious caveat with the detection is the lack of flat base- 
line in the spectrum above « 300 km s _1 due to the lost chan- 
nels. However, the spectra obtained in each individual polarisation 
are fitted by very similar Gaussians. Furthermore, if instead we use 
3C 286 for the bandpass calibration, which was observed at the start 
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Velocity relative to z=0.6561 [km s _1 ] 

Figure 1. Top: The bandpass calibrated spectrum of 3C 286. Middle: The 
21-cm spectrum of 1622+238. The Gaussian fit gives a peak line depth of 
30.8 mJy (giving r = 0.0088, against a 3.55 Jy continuum flux density), a 
centroid at -47±5 km s _1 and a FWHM of 235 km s -1 for the line. The 
centroid gives a frequency of 857.815 MHz, which corresponds to a redshift 
of 

^abs — 0.65584 i 0.00003, cf. the values determined for the Lyman-a 
absorption: 0.6561 (Rao & Turnshek 2000) and 0.6564 (Chen et al. 1998). 
The r.m.s. noise is 4.7 mJy per each 10 km s" 1 channel. Bottom: The ab- 
sorption of Mgl (unsaturated) and Mgll (saturated) at z = 0.656 towards 
1622+238 (from ESO service-mode program 69.A-0371). 

of the run, rather than 3C 48, which was observed at the end, the 
data give very similar profiles (again with the 1 Jy dip in the flux 
below 856.6 MHz), which are fitted by the same Gaussian parame- 
ters as before. As would be expected from a real feature towards an 
unresolved (< 75" x 42") background source, the line is present 
in the centre and surrounding pixels, but disappears with increas- 
ing proximity from the centre pixel. We also produced a cube in 
which 1622+238 is offset from the phase centre: While the absorp- 
tion was still apparent towards 1622+238, no feature was seen at 
the phase centre. As a further check, we used 3C 48 to bandpass 
calibrate 3C 286, and, as seen from Fig. 1 (top panel), there is no 
overall ripple in the spectrum. 

3 RESULTS AND DISCUSSION 
3.1 Observational results 

At t = 0.0088, the peak optical depth is at the lower end of 
the DLAs detected in 21-cm absorption (table 1 of Curran et al. 
2005), although the large line-width gives a velocity integrated op- 
tical depth of J rdv — 2.2 ± 0.2 km s~\ which takes it to joint 



Table 1. The DLAs detected in 21-cm absorption listed in order of increasing 
spin temperature/covering factor ratio. 2 a b s is the redshift of the DLA, J rdv is 
the velocity integrated optical depth of 21-cm absorption [km s -1 ], TVhi is the 
total neutral hydrogen column density [cm~ 2 ] of the absorber (see Curran et al. 
2002), which is unknown for the DLAs towards 0248+430 and 2351+456. The 
optical identification of the absorber (ID) is given: D-dwarf, L-LSB, S-spiral, 
U-unknown. T s /f is the inferred spin temperature/covering factor ratio [K]. The 
final column gives the 21-cm absorption reference. See Curran et al. (2005) for 
details. 
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5th place in the line strength stakes (Table 1). In the optically thin 
regime (r < 0.3), the total atomic hydrogen column density, A^hi 
[cm -2 ], of the absorbing gas in a homogeneous cloud is related to 
the velocity integrated optical depth of the 21-cm line via 

AThi = 1-823 x 10 18 y J t dv , (1) 

where T s [K] is the spin temperature of the gas, / is the cover- 
ing factor of the flux by the absorber and J rdv is the velocity 
integrated optical depth of the 21-cm absorption. If the 21-cm and 
Lyman-Q absorption arise in the same cloud complexes, the total 
neutral hydrogen column density of 2.3 x 10 20 cm -2 (Rao & Turn- 
shek 2000), gives T s /f ~ 60 K, the lowest value yet detected in 
a redshifted source (Table 1). Since both the spin temperature and 
covering factor are generally unknown, the degeneracy in the ra- 
tio is best left intact (Curran et al. 2005; Curran & Webb 2006), 
although for / ^ 1 (by definition) this means that T s ^ 60 K. 

The 21-cm absorption profile in this absorber is much wider 
than that found in other DLAs ($: 53 kms -1 , see table 1 of 
Curran et al. 2005) and more like the typical Hi emission pro- 
file of a galaxy (e.g. Koribalski, et al. 2004). We note also, in 
the other DLAs the 21-cm profile-width is always narrower than 
the metal-ion lines (typically FWHlvLicm ~ 0.2 AVivigii, Curran 
et al. 2007), but in this case FWHM21™ ~ 2 AVu s U, which sug- 
gests very different kinematics from those of the typical DLA. In 
fact, for other DLAs the spin temperature is often compared with 
the kinetic temperature (e.g. Lane & Briggs 2001), which should 
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be equivalent for single cloud in thermodynamic equilibrium. Ap- 
plying T k « 22 x FWHM 2 to 1622+238 gives the unphysical 
Tk ~ 10 6 K, confirming that kinematical motions are chiefly re- 
sponsible for the wide profile. 



3.2 The nature of the z = 0.656 absorber towards 
1622+238 

3.2.1 Possible absorber hosts 

In Fig. 2 we show the HST WFPC2 (F702W filter) optical image 
of the field towards 3C 336 (z om = 0.927) overlaid with the radio 
maps. Although there is compact 8.4 GHz emission from the QSO 
itself, most of the radio emission is from the extended lobes, which 
span an enormous 28" (Bridle et al. 1994; Cheung et al. 2005). 
The DLA towards this QSO is somewhat of a mystery, since the 
only galaxy in the HST field with a redshift of z a b s = 0.656 must 
have a projected disk impact parameter of « 65 hT 1 kpc (Steidel 
et al. 1997). Since no other galaxy within this impact parameter 
has ever been found, Steidel et al. (1997) suggest, but stipulate as 
unlikely, that the absorption may be due to an under-luminous late- 
type spiral, in which an invisible disk spans at least all of the way 
from the identified z = 0.656 galaxy to the QSO sight-line. 

As an alternative, the authors suggest that there may be other 
galaxies within this impact parameter, although as these are unde- 
tected, they would have to be extremely faint (L < 0.05L^). An- 
other possibility is offered by the galaxy of unknown redshift close 
to the optical hotspot, or even the hotspot itself, both of which are 
located close to a radio lobe sight-line (Cheung et al. 2005). If these 
(and perhaps some other low luminosity galaxies), were part of a 
group, which included the z — 0.656 galaxy, the redshift criterion 
would be satisfied, while perhaps providing a large H I reservoir lo- 
cated closer to the QSO sight-line. Furthermore, the H I could be 
largely dispersed over the region through the mutual tidal interac- 
tions of these galaxies (Steidel et al. 1997). 



3.2.2 Implications of the 21-cm absorption 

From Fig. 2, we see that, although the identified z = 0.656 galaxy 
is located close to the sight-line of a radio lobe, it remains suf- 
ficiently remote to be covering a negligible amount of flux 1 . Al- 
though we cannot determine the location of the absorbing region 
within the radio-bright region, we may obtain an estimate of its 
extent: The circular region centered on the z — 0.656 galaxy rep- 
resents the physical area which encompasses 0.88% of the total 858 
MHz flux. That is, this is the smallest region that can give the ob- 
served optical depth. So although the visible feature identified as 
the z = 0.656 galaxy does not intercept much radio flux, a cold, 
invisible disk of Z 20 kpc (3.3") radius could cover significant 
radio emission (as shown in the figure), although, as per the reser- 
vations of Steidel et al. (1997), it may not necessarily be centered 
upon this galaxy. In fact, the 21-cm and Lyman-a absorption may 
not be coincident with one another, in light of the relative profile 
widths and offset in the profiles (Fig. 1). 

For a zero intrinsic velocity offset, one can estimate the vari- 
ability of the constant x = a 2 g p fi, where a is the fine structure 



1 Although we would expect the 1420 MHz emission in the rest frame 
of the galaxy (858 MHz) to be further extended than the 21-cm emission 
measured in our rest frame. 



constant, g p the proton (/-factor and /i = m c /m p is the electron- 
to-proton mass ratio. In such a case, the relative change in the value 
of x between redshift z and a terrestrial (z — 0) value is given by 
(e.g. Tubbs & Wolfe 1980) 

Az _ X z - Xp _ Z V y - Z 2 l 
X X 1 + Z21 

where 221 and ztjv represent the observed absorption redshifts 
for 21-cm and a UV heavy-element transition, respectively. Us- 
ing the Gaussian fit centroid, 221 = 0.65584 ± 0.00003, and the 
stongest Mg I component, z uv = 0.656062 ± 0.000003, we ob- 
tain Ax/x = (1.34±0.18) x 10~ 4 .This is an order of magnitude 
larger than any of the eleven other values which have been deter- 
mined (see figure 1 of Tzanavaris et al. 2007), suggesting that the 
velocity offset between the 21-cm and UV absorption is not primar- 
ily due to the values of the fundamental constants at z — 0.656. 

Therefore the observed w 50 km s _1 velocity difference, in 
conjunction with the atypical 21-cm line-width, suggests a physical 
offset between the cold, 21-cm absorbing, and warm, UV absorb- 
ing, gas. However, although not completely coincident, hydrogen 
is observed in both regimes, within the same field (possibly remote 
from the QSO sight-line), at (nearly) the same redshift, and so we 
find it unlikely that the radio and UV absorption are not related 
in some way. The width of the 21-cm absorption in comparison to 
the optical lines, and that of other DLAs, may be the result of an 
atypical situation where we have an extended absorber occulting 
an extended radio source. However, all of the known DLAs associ- 
ated with spiral galaxies also occult large radio sources (figure 4 of 
Curran et al. 2005). Furthermore, although located much closer to 
the spirals than the other DLAs in terms of T s / f, due to the wide 
profile, the 21-cm absorption towards 1622+238 may be somewhat 
of an anomaly (Fig. 3) 2 . 

If the system is bound and the wide 21-cm profile is purely 
the result of rotational broadening (i.e. not an outflow), the total 
line-width of ±280 kms -1 , leads to a dynamical mass estimate 
of M dyn = 1.8 x 10 1() r/sin 2 i M , where r [kpc] is the edge 
of the cold H I distribution [or the 858 MHz emission, whichever 
comes first]. In the case where this gas is located close to the QSO 
sight-line, so that the identified z = 0.656 galaxy is actually near 
the edge of the absorbing disk or group 3 , i.e. r m 75 kpc (Chen 
et al. 1998), the dynamical mass becomes M dyn ~ 1 x 10 12 / sin 2 i 
Mg. If the gas is due to a large spiral, centered on the z = 0.656 
galaxy 4 , by assuming circular motions and the same inclination as 
the visible "nucleus" (55°, Chen et al. 1998), the deprojected ve- 
locities and impact parameter give Md yn ~ 2 x 10 12 M Q (for 
r ~ 75 kpc). Although large, the derived masses are close to 
the values derived for the two most local spirals, the Milky Way 
(Mdyn ~ 2 x 10 12 M Q , Wilkinson & Evans 1999) and M31 
(M dyn ~ 1 x 10 12 M , Evans et al. 2000), and certainly rea- 
sonable for a group of smaller galaxies. 

Combining this dynamical mass with the luminosity of L < 
0.25 L* (Steidel et al. 1997) 5 , gives a mass-to-light ratio of 



2 Note that if FWHM 2 i cm ~ 0.2 AV MgII , as per the other DLAs, T B /f 
would be Ri 600 K, placing this bang in the middle of the spirals in Fig. 3. 

3 Assuming i is the inclination of the most rapidly orbiting galaxies and 
neglecting their individual rotations. 

4 Which may be indicated by the single Gaussian fit, although signal-to- 
noise ratio of the spectrum (Fig. 1) is far too low to resolve possible contri- 
butions from smaller discrete galaxies. 

5 Note that, from local studies, the vast majority (87%) of DLAs are be- 
lieved to be due to galaxies of luminosities < L* (Zwaan et al. 2005). 
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Figure 2. The central part of the HST image of Steidel et al. (1997) overlaid with the 8.4 GHz [inner/red contours] and 1 .4 GHz [middle/green] radio emission 
(Cheung et al. 2005). We also show the FIRST 1.4 GHz radio map [outer lobe/blue]. We do not show the 858 MHz continuum emission from our observations 
as the quasar is unresolved by the 75" X 42" WSRT beam. The contours scale from 10% of the peak emission, and for the sake of clarity, shown only up to 
50% in the 1.4 GHz maps. The only z = 0.656 galaxy found by Steidel et al. (1997) is highlighted. The circular region, centered at this galaxy, shows the 
size of a region which contains 0.88% of the total 858 MHz flux density (not necessarily located here, see text). The linear distance indicated is calculated for 
a redshift of z = 0.656 using Ho = 75kms~ 1 Mpc -1 , f2 ma ttcr = 0.27and Oa = 0.73 (used throughout the paper). North is up and east to the left. The 
HST image and radio maps are courtesy of Teddy Cheung and the Very Large Array's FIRST (Faint Images of the Radio Sky at Twenty Centimetres) survey. 



M dyn /L > 400 M q /Lq. Such ratios (~ 1000 M Q /L ), define 
the hypothesised low density galaxies (LDGs), which can range in 
mass from 10 9 to 10 12 M Q (Jimenez et al. 1997). In an LDG a low 
density disk resides in a dark halo, with the surface density be- 
low that required by the Toomre criterion for instabilities to form 
(^ 10 20 cm -2 , Jimenez et al. 1997; Verde et al. 2002), and so there 
is no ongoing star formation: From deep HST Ha imaging, Bouche 
et al. (2001) determine a star formation rate of < 0.15 Mq yr -1 
kpc - in the 3C 336 field. Another possible example of an LDG 
may already have been observed as the "dark hydrogen cloud", 
VIRGOHI 21 in the Virgo Cluster (Minchin et al. 2005). This has a 
mass-to-light ratio of M^ yn /LB > 500 Mq/L© and spans at least 
16 kpc in extent. If it is a galaxy, this has a mass of Afdyn > 10 



Mq, although, as may be the case towards 1622+238, the large 
amount of gas could be due to tidal debris. Where ever the H I (Ly- 
a & 21-cm) absorption arises towards 1622+238, it is clear that 
this would not be the first time a large amount of unidentified neu- 
tral hydrogen has been detected. 



4 SUMMARY 

We have detected 21-cm absorption in the z a b s = 0.656 damped 
Lyman-Q absorber towards 1622+238 (3C336). If the absorption 
arises in the same gas as the Lyman-a absorption, the spin tem- 
perature of the gas is T s ^ 60 K, the lowest yet found in a 
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Figure 3. The spin temperature/covering factor ratio versus the absorption 
redshift for the DLAs searched for in 21 -cm absorption. The symbols repre- 
sent the 21 -cm detections and the shapes represent the type of galaxy with 
which the DLA is associated: circle-unknown type, star-spiral, square- 
dwarf, triangle-LSB. The arrows show the lower limits and all of these bar 
one (0454+039 at 2 a b s = 0.8596) have unknown host identifications. The 
unfilled symbols show the new detections since Curran et al. (2005), from 
which the figure is adapted: Kanekar et al. (2006, 2007) and 1622+238 (un- 
filled star), on the assumption that the 21 -cm absorption arises in the same 
gas as the Ly-a absorption. 

DLA. The low spin temperature is a consequence of the extremely 
wide profile, which, most atypically, is wider than the Mg II line. 
This could be due to a large absorber occulting what we know to 
be large radio source. The 21 -cm is also offset from the unsatu- 
rated Mgl absorption by ~ —50 kms^ 1 , which corresponds to 
Ax/x »s 1.3 x 10~ 4 (x = a 2 g p jj,), which is an order of magni- 
tude larger than that expected from other DLAs, suggesting that the 
offset is not dominated by different values of the fundamental con- 
stants at z = 0.656. It may indicate, instead, that the radio and UV 
absorption are unrelated, although a lower Lyman-a column den- 
sity would give gas significantly colder than 60 K. Furthermore, 
it seems somewhat strange that the hydrogen absorbers detected 
in two different regimes at the same redshift would be unrelated, 
especially considering that none of these seems to arise directly 
along the QSO sight-line. Regarding this, our results cannot distin- 
guish from where the hydrogen absorption arises, whether this be 
the disk of a large underluminous spiral, an extremely faint galaxy 
lying close to the QSO sight-line or a possible tidal interaction, be- 
tween other faint galaxies, which is dragging significant amounts 
of gas across the sight-line. 
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